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ABSTRACT 

Recent Chandra observations of interacting and starburst galaxies have led us to investigate the 
apparent correlation between the positions of young star clusters and Chandra point sources. Assumed 
to be X-ray binaries (XRBs), these point sources do not seem to coincide with the massive (~10^ Mq), 
young (1 — 50Myr) stellar clusters that can easily form systems capable of such emission. We use 
a sophisticated binary evolution and population synthesis code (StarTrack) and a simplified cluster 
model to track both the X-ray luminosity and position of XRBs as a function of time. These binaries 
are born within the cluster potential with self-consistent positions and velocities and we show that 
a large fraction (~70%) can be ejected from the parent due to supernova explosions and associated 
systemic velocities. For brighter sources and cluster masses below ^IO^Mq, we find that the average 
number of bright XRBs per cluster remains near or below unity, consistent with current observations. 

Subject headings: galaxies: star clusters - methods: statistical - X-rays: binaries 



1. INTRODUCTION 

One of the interesting developments made possible 
with Chandra observations is the abundance of extra- 
galactic point sources. The large observed samples 
now available promise to improve our understanding of 
their formation and interaction with their host stellar 
environments. In galaxies with young stellar popula- 
tions, it is generally thought that brig ht point X-ray 
sourc es are young X -ray Bina ries fXE Bsl (iKilgard et al I 
l2(ffll IFabbiano et a l. '2001* iFabbiano fc White 1 1200?^ 
associated with the large amount of ongoing star for- 
mation. This interpretation is mainly based on the 
measured X-ray luminosity and the spectral and tem- 
poral variability characteristics of the point sources. 
For a detailed re view of the formation and evolu- 
tion of XRBs see iLewin fc van der Klis I l(2005(l . and 
for a recent discussion of th e te mpor al properties see 
iSipior. Eracleous. fc Sieurdsson I <j2003|l . Optical and in- 
frared observations, most prominently with the Hubble 
Space Telescope, reveal massive, young clusters, often 
referred to as super star clus ters. T hey range in mass 
from^~10^M(T, to ~10^Mp (jSmith fc Gallagher 200ll 
Harris et al. 1 120 011) . and in age from just ^IMyr to 
~50Myr, with the majority at ~ 10 — 20Myr (mostly due 
to photometric selection). These clusters are thought to 
be young analogs of globular clusters and may be respon- 
sible for most of the massive star s in the field of their host 
galaxies (jTremonti et al. ir200H) . Thus, one may expect 
a concetrati on of XRBs in or nea r these clusters. 

Recently, IKaaret et al. I (j2004() : hereafter K04 have 
studied 3 starburst galaxies (M82, NGC1569, and 
NGC5253) each containing a significant number of young 
star clusters and point X-ray sources. Indeed, they do 
find a statistically significant relationship between the 
two types of objects: XRBs are preferentially found 
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within distances of ~30— lOOpc from their nearest clus- 
ter, but there is a clear lack of XRBs found coincident 
with the clusters. There are obvious observational biases, 
the most important being that the true parent cluster 
is unknown; the distances quoted are only those to the 
nearest cluster, not necessarily the parent cluster. Still, 
the XRB spatial distribution relative to the clusters and 
their association seems significant and characteristic of a 
non-random sample distribution (see § 3 of K04 for the 
relevant statistical analysis). It is worth noting that sim- 
ilar res ults have been fou nd in observations of the An- 
tennae l|Zezas et al. II200I . 

In this Letter, our goal is to model in a self-consistent 
manner the population of binaries in the cluster potential 
(assumed static for simplicity). We track the kinematic 
evolution of compact object binaries in the absence of dy- 
namical interactions and we follow their X-ray luminosity 
(Lx) evolution. We focus on two specific, testable points 
of comparison between the published observations and 
our calculations: the average number of XRBs per clus- 
ter, and the median distance of XRBs from their parent 
cluster (or nearest cluster in the observations). For point 
sources brighter than ~10'^^ erg s~^ , K04 have shown that 
the median distance from a cluster is ~100pc, with an 
average of <1 XRB per cluster (See Table 1 in K04 for the 
specifics regarding each of the three galaxies considered). 
Here we show that these two quantities can be calculated 
theoretically and the results appear consistent with the 
observations for the range of cluster masses and ages rel- 
evant to the clusters observed, when considering only the 
supernova kicks imparted to XRBs and their motion in 
cluster potentials. Guided by the Lx sensitivity limits 
of the observations, we focus our analysis on XRBs with 
Lx >5xl0^^ergss~"'^, although the formation and evo- 
lution of XRBs to lower Lx ranges is included in our 
models. 

In § 2 we describe the model methods used and how 
they are applied in our simulations. In § 3 we describe 
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our main results and compare them to the observations 
presented by K04. We discuss our conclusions in § 4. 

2. THEORETICAL MODELING 

To generate the necessary stellar populations for the 
modeled clusters, we use the population synthe s is pro - 
gram StarTrack (developed bv iBelczvnski et al~l l)2002t) : 
Belczynski et al. 2004, to be submitted). We generate 
and evolve a population of binaries under a given set of 
conditions, such as the initial mass function (IMF), su- 
pernova kick distribution, common envelope efficiency, 
etc. With the resultant evolutionary parameters of the 
binaries at the time of the compact-object (neutron star 
or black hole) formation, we place them in a cluster po- 
tential and track their motion and X-ray luminosity as 
a function of time. In so doing, we ultimately generate 
a complete evolutionary picture of the X-ray binaries in 
association with their parent cluster. As noted already, 
we do not account for any stellar interactions in these 
young clusters, as our goal is to examine whether su- 
pernova kicks alone can account for the observed spatial 
distribution of XRBs relative to their parent clusters. 

StarTrack is a sophisticated Monte Carlo population 
synthesis code that has been recently updated to care- 
fully account for binary mass-transfer phases and Lx 
calculation. We account for various phases of mass and 
angular momentum losses, and have implemented an in- 
tegrated tidal evolution method that is calibrated against 
observations of Galactic high-mass XRBs and of circu- 
larization in open clusters. Some key features for this 
investigation include: (i) the determination of the post- 
core-collapse systemic velocity for compact object bina- 
ries, and (ii) the detailed calculation of the mass transfer 
rate between binary components, calibrated against cal- 
culations with a stellar evolution code. Systemic veloci- 
ties are a key to the proper determination of the orbital 
trajectory, which is one of the primary concerns for this 
work. Also, given the sensitivity of the observations to 
Lx, the mass transfer rate becomes a critical factor in 
determining whether or not a given XRB is relevant to 
the K04 observations at any point in its lifetime. For the 
Lx determination we further apply a bolometric correc- 
tion to the theoretical value, to account for Chandrds 
sensitivity band and the typical XRB spectra. This bolo- 
metric correction is dependent on the system parameters 
(neutron star or black hole accretor, wind accretion or 
Roche-lobe overflow) and assumptions about the typi- 
cal spectra of different sources, derived empirically from 
Galactic observations (Maccarone 2003, private commu- 
nication; iPortegies Zwart. Dewi. fc Maccaronel 1)2004(1 '). 
Specifically for wind accretors we adopt 0.15 and 0.7, for 
disk persistent sources and transient sources at outburst 
we adopt 0.5 and 0.7 (for neutron stars and black holes, 
respectively) . 

It is well known that population synthesis calcula- 
tions require a significant number of input parameters. 
Since our main goal in this Letter is a proof-of-principle 
study, it is not necessary to fully explore the param- 
eter space. Instead we choose to consider a reference 
model with parameter assumptions that are considered 
ty pical for binary e v olutio n calculations (see model A 
in IBelczvnski et al. I l|2002fl 'l. We also consider a small 
set of other models where we vary the initial mass func- 
tion of binary primaries, as this most significantly affects 



TABLE 1 

Parameters for model runs 



Model # of Mass IMF Ri/j 
MC runs (Mq) index (pc) 



A 


1000 


SxlO"* 


2.35 


10 


B 


1000 


5x10'' 


2.7 


10 


C 


100 


5x10^ 


2.35 


10 


D 


100 


5x10^ 


2.7 


10 


E 


7 


5x10^ 


2.35 


10 


J 


500 


5x104 


2.35 


1 



the relative contribution of XRBs with neutron stars and 
black holes that acquire systematically different systemic 
velocities. 

The systemic velocities are determined by the na- 
tal kicks imparted during the formation of the 
compact object. Neutron star natal kick magni- 
tudes are drawn from the dist ribution derived from 
lArzoumanian. Chernoff. fc Corde s ( 2002) , based on cur- 
rent pulsar kinematics. It consists of two Maxwellians 
with cr's 90 and 500, with a relative weight of 2:3, re- 
spectively. Black hole natal kicks are linearlly scaled 
from Neutron star kicks based on their mass relative to 
a typical Neutron star mass of 1.44 Mq. 

For the orbital evolution we use the calculated post- 
core-collapse systemic velocities of XRB progenitors and 
combine them consistently with a Plummer model for the 
gravitational potential of the model cluster with a given 
half-mass radius, Ri/2- To determine the spatial distri- 
bution of the XRB progenitors right after compact-object 
formation (which generally occurs immediately prior to 
the X-ray phase) , we assume that the number density of 
stars is proportional to the mass density. Initial systemic 
velocities consi stent with the P l umme r potential are also 
generated (see lAarseth et all 1)19741) V We then apply 
the calculated post-core-collapse systemic velocities ran- 
domly oriented with respect to the initial cluster veloc- 
ities. We follow the motion of the binary as a function 
of time and correlate position with the X-ray luminosity 
evolution calculated with the binary evolution code. 

The number of binaries modeled for a given cluster 
is directly proportional to the mass of that cluster once 
the IMF index, mass-ratio-distribution paramters, and 
binary fraction are chosen. We adopt a flat mass-ratio 
distribution and a binary fraction equal to unity in order 
to represent an upper limit on the number of binaries, 
and therefore on the total number of X-ray sources. We 
consider cluster masses in the range 10^ M© — 10^ M©, 
with IMF indices of 2.35 and 2.7. The specific parameters 
of our simulations are shown in Table 1. 

It is interesting to note that although orbits are calcu- 
lated typically for few to several hundred Myr, binaries 
are X-ray sources for only a small part of their orbit and 
they are bright (i.e., Lx >5x 10'^^ ergs~^) X-ray sources 
for an even smaller part. Each system is evolved in- 
dividually in a static cluster potential, and thus no in- 
teractions or cluster evolution is allowed in the present 
analysis since our goal is to examine the effect of the 
supernova kicks. This may not be a well justified as- 
sumption in general. Nevertheless, the clusters relevant 
to our study are very young (few Myrs to ~ 10 — 20 Myr), 
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so significant cluster evolution is not expected, except for 
possibly in the most massive and most compact (small 
half-mass radius) clusters. 

We note that statistical effects play a significant role, 
especially in the low mass ('^lO^'M©) clusters. Typi- 
cally, no more than one XRB is bright enough to be 
seen in these clusters and the position can vary signif- 
icantly across the cluster for each different simulation. 
Therefore, we consider a large number of Monte Carlo 
realizations for each parameter set (cluster mass, half- 
mass radius, and IMF index). The lower mass clusters 
have the smallest number of initial binaries and hence 
require the most realizations. We chose the number of 
realizations for each cluster mass so that our results av- 
eraged over the many realizations remained unaffected 
at the 5% level. 

3. RESULTS & DISCUSSION 

Our calculations have yielded a wealth of results. 
Knowing both the trajectory and X-ray luminosity of 
such a wide range of objects can help us understand 
XRB formation in young clusters both statistically and 
in a system-by-system sense. Here, we focus on the sta- 
tistical results where the clusters are investigated as a 
grand average; this is more appropriate when comparing 
with large populations of clusters (as analyzed in K04). 
Specifically, we focus on the statistical averages of the 
two quantites quoted observationally in K04: the me- 
dian distance of XRBs from the nearest (parent in the 
models) cluster, and the mean number of bright XRBs 
per cluster (within 1000 pc). 

In Figure 1 (top) we plot the model average number 
of XRBs per cluster as a function of distance from the 
parent clusters, each of 5x10'* Mq and for a variety of 
cluster ages. These ages are within the range of estimates 
for the observed clusters. To take into account the un- 
certainties in the age estimates (typically a few Myr), we 
use an age-snapshot method, based on which we deter- 
mine the average number of XBRs as a function of radius 
for a specific "instant" in time, and then average these 
results over each consecutive "instant" within the cluster 
age estimate and its error. 

It is evident from Figure 1 that the XRB spatial distri- 
butions have a dramatic time dependence. For "young" 
clusters, the average XRB number per cluster rises to 
a maximum rapidly and very few XRBs are found at 
large distances. This is primarily because even the un- 
bound XRBs have not had enough time to move away 
from their parent clusters. The median systemic veloc- 
ity of the XRB systems is ^2— 3pcMyr~*, which limits 
the distance any XRB can reach. For older clusters, the 
average XRB number exhibits a fairly slow increase with 
distance, up to 2kpc and sometimes beyond. This can 
potentially create a pollution effect and lead to difficulty 
in identifying the true parent cluster in observations. 

It is also evident that at certain ages XRBs are 
distinctly more numerous than at others. For exam- 
ple, in Figure 1, the SxIO^Mq (top) clusters with 
Lx >5xl0^^ergs~^ show more XRBs at 10 Myr than 
at any other time in the clusters' evolution. We also find 
that this peak age is dependent on the Lx cut-off. Fully 
exploring these dependencies could allow us to derive 
general conclusions about XRB populations dependent 
only on the average ages of the young cluster popula- 




Distance from Cluster [pc] 

Fig. 1. — Average number of XRBs seen within a given distance 
from its parent cluster for specific cluster ages listed shown. See 
Table 1 for the specifications of the model parameters. All data 
shown here is calculated with an Lx cutoff of 5x10^^ ergs~^. 



tion. 

In Figure 1 (middle and bottom) we present our results 
for clusters of 5x10^ and SxlO^M©, respectively. Note 
that the behavior is similar for all masses, except that 
the average number of XRBs at a given radius scales 
with the mass of the cluster almost lineraly. This is due 
to the direct relationship between the number of binary 
systems modeled and the cluster mass. 

We calculate the median distance and mean number 
of XRBs with Lx >5xl035ergs-^ within 1000 pc (Ta- 
ble 2), in order to compare appropriately with K04. 

Mean number of XRBs per cluster: We find the the- 
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_ TABLE 2 

Mean XRB Number {Nxrb) and Median^ XRB distance from cluster center (Rmedian) 







5 Myr 


10 Myr 


15 Myr 


20 Myr 


25 Myr 


50 Myr 


Model 


X RB ^median 


XRB 


^median 




^median 


J^XRB 


'^median 


J^XRB 


"^median 


Nxrb 


-^median 






(pc) 




(pc) 




(pc) 




(pc) 




(pc) 




(pc) 


A 


0.19 


10.5 


0.23 


17.5 


0.18 


44.5 


0.18 


54.5 


0.19 


58.5 


0.04 


146.5 


B 


0.14 


10.5 


0.18 


17.5 


0.12 


46.5 


0.12 


67.5 


0.16 


76.5 


0.04 


104.5 


C 


2.00 


9.5 


2.55 


15.5 


1.72 


27.5 


1.87 


46.5 


2.26 


51.5 


0.41 


80.5 


E 


20.5 


16.5 


24.3 


17.5 


17.3 


18.5 


18.1 


24.5 


20.4 


26.5 


6.00 


29.5 


J 


0.18 


1.5 


0.25 


2.5 


0.18 


13.5 


0.16 


31.5 


0.21 


40.5 


0.04 


54.5 



^Only XRBs within 1000 pc are used to calculate the values listed here in order to compare with K04 



oretical mean XRB number per cluster to vary signif- 
icantly from ~ 0.1 to ~ 10, depending on the cluster 
mass. Therefore, it is possible to reproduce the results 
in K04 by taking contributions from a number of clusters 
of different masses. Two of the three galaxies discussed 
in K04 (M82 and NGC 5253) have a mean number of ob- 
served XRBs of ~ 1 per cluster, while NGC 1569 seems 
to have a very small number of XRBs (only ~ 0.25 per 
cluster). This difference would point towards NGC 1569 
having, on average, smaller-mass clusters, even though 
outliers at high masses can still exist. A difficulty in the 
comparison arises because the properties of the clusters 
in these galaxies are difficult to determine orbservation- 
ally. Those with measured masses are skewed to higher 
masses (>1 x lO'^ Mq) and younger ages (<15Myr) sim- 
ply because they are selected photometrically (Gallagher 
2004, private communication). Therefore, developing a 
proper theoretical cluster distribution for comparison is 
rather challenging without further observational studies 
of the cluster populations. 

Median distance of XRBs from the cluster: Our results 
(Table 2) indicate a strong dependence of the median 
XRB distance on the age and a moderate dependence 
on the cluster mass. For clusters with a half-mass ra- 
dius of 10 pc and masses ^ 5 x 10^ M©, median distances 
reach values of 30 — 100 pc (similar to those observed) 
at ages of 15 Myr and older. Only very massive clusters 
of ^ 5 X 10^ Mq reach such distances later at ~ 50 Myr. 
These ages and moderate masses are consistent with the 
current observational estimates, although massive and 
older clusters are also present in the photometrically se- 
lected clusters in K04 (Gallagher 2004, private commu- 
nication) . 

It should also be noted that, for the highest cluster 
mass we consider (5x10® M0), even the oldest clusters 
seem to show more binaries than what is observed. This 
clearly implies that starbursts are not dominated by such 
massive clusters, and this is not surprising. However, 
these more massive clusters may also be affected by dy- 
namical cluster evolution and stellar interactions leading 
more binary disruptions and ejections. Thus we would 
expect the average number of XRBs per cluster to de- 
crease at all ages and distances. 

It should be noted that we have assumed a binary frac- 
tion of unity, and therefore the mean XRB numbers could 
be overestimated. This is true also because projection ef- 
fects have not been taken into account, and our numbers 
represent the radial distance the XRBs have traveled. 
Also, we note that changes in the power law IMF index 



of the cluster produce noticeable, but largely insignificant 
changes in the cluster profiles. For example, changing the 
IMF index from 2.35 to 2.7 decreases the average number 
of binaries at or about the 10% level for each timestep. 
This effect may become more important, especially for 
very steep IMFs, su c h as t hose for clusters proposed by 
iKroupa fc Weidnerl l)2003|) where the index can go as 
high as 3.2. And last, changes in the half-mass radius 
of the cluster dramatically change the median XRB dis- 
tance for a given mass. Very small values (model J in 
Table 2) tend to limit XRBs ejection, as the potential is 
deeper. In these tight clusters, it is likely that dynamics 
will play a non-negligible role, depending on their age. 

4. CONCLUSIONS 

With detailed population simulations of XRBs and 
a simple treatment of gravitational potentials of young 
clusters we have shown that the significantly low XRB 
numbers per cluster observed in starbursts can be ex- 
plained as being largely due to supernova kicks im- 
parted to XRBs at compact-object formation that 
lead to XRB ejection from the cluster potential, as 
heuristically suggested by observational stu dies (K04; 
iPortegies Zwart. Dewi. fc Maccaroiiel ()2004D ). Derived 
XRB median distances are also consistent with current 
estimates of cluster masses and ages, although a more 
direct comparison requires more detailed observational 
constraints of the cluster properties. 

This work opens many possible avenues in which to 
continue this study, some of which include an in depth 
look into the systematics generated by our stellar evolu- 
tion code, such as how our results change with a broader 
range of masses and IMFs, as well as additional stel- 
lar evolution parameters such as the common envelope 
efficiency. We also intend to look at the detailed popula- 
tions created, and search for specific correlations between 
types of XRBs, their ages, and positions in the clusters. 
Still further, we have largely ignored the low luminosty 
XRBs in this analysis. This population may indeed be 
detectable, if present in large enough quantities, as dif- 
fuse emission. And lastly, it is possible that for the more 
massive, compact, and older clusters, dynamics play a 
non-negligible role in the XRB evolution. We hope to 
extend our modeling to include dynamical considerations 
such as this in the near future. 

We are grateful to J. Gallagher, T. Maccarone, and 
P. Kaaret for useful discussions. The work is partially 
a Packard Foundation fellowship and a NASA Chandra 
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